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Abstract. A poorly constrained parameter in the Standard Model of Cosmology is the lepton 
asymmetry I = If = ~^ n^^)/s. Each flavour asymmetry If with / = e, r is the 

sum of the net particle density of the charged leptons Uf and their corresponding neutrinos, 
normalized with the entropy density s. Constraints on Z/ from BBN and CMB allow for lepton 
flavour asymmetries orders of magnitudes larger than the baryon asymmetry b ~ 10"^". In this 
article we show how large lepton (flavour) asymmetries If < 0(0.1) influence the early universe, 
in particular the freeze out of WIMPs and the cosmic QCD transition. 



1. Introduction 

The recently reported higher primordial abundance of ^He |lj and the larger number of effective 
number of neutrinos by observations of the cosmic microwave background [21 |3] might be a 
motivation for large neutrino chemical potentials or equivalently large neutrino asymmetries. 
Charge conservation and neutrality ensures that the asymmetry in the charged leptons is of the 
same order as 6 = nb/s{T) ~ 0(10^^'^) [2]. But due to the difficulties in detecting neutrinos 
one can only obtain constraints on the lepton (flavour) asymmetry from indirect observations. 
A large If would affect big bang nucleosythesis (BBN) in two ways. It changes the primordial 
■*He abundance [Ij and the expansion of the universe [5]. We show in this article some results 
on how large lepton (flavour) asymmetries influence the early universe between the electroweak 
phase transition and the onset of neutrino oscillations. For a more detailed discussion see [6] 
and 17|. 

2. Pre-BBN 

Before BBN, at Tqsc — few MeV, neutrino flavour starts to oscillate and equilibrates any initial 
flavour asymmetries before the onset of BBN [8j. For temperatures larger then Tqsc lepton 
flavour number is conserved. Within the Standard Model, no violation of baryon and lepton 
number are possible after the electroweak phase transition at Tew ~ 200 GeV [9|. A lepton 
flavour asymmetry If remains constant in the early universe between Tew and Tosc- Here we 
assume // < 0(0.1) HE]. 

All particle interactions between these two events are in statistical equilibrium and can be 
described by the set of 5 equations: 






b s(T) = ^ BiUb with Bi = baryon number of species i (2) 

i 

q s{T) = = ^ QiUq with Qi = electric charge of species i. (3) 

i 

The net particle density of a species i with mass and chemical potential /Zj is given by 
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^ , E^E^ - rrii^ p \dE (4) 

with the internal degrees of freedom of a species gi. Anti-particles i have the same chemical 
potential as the corresponding particle, but inverted sign: /ij = — With the observed baryon 
asymmetry and the charge neutrality of the universe we can solve equations ([T]) to ([3]) for different 
1 1 for all Standard Model particles in statistical equilibrium, leaving \Xi as free parameters. 

With each conserved quantum number one can associate a chemical potential which can be 
calculated from the free energy of the particle plasma. Apart from the three lepton chemical 
potentials and the charge chemical potential we can also calculate the baryon chemical potential 
/Ufo. The latter one is for cosmologist and particle physicist of greater interest to describe the 
QCD transition, the confinement of quarks to hadrons. It is still an open question, if the 
cosmic QCD transition is just a rapid change or a real phase transition. One would like to have 
deeper knowledge of the order of the transition, since it sets the initial conditions for BBN and 
the different orders would lead to different observable consequences, like a modification of the 
primordial garavitational background or quark nuggets [11] . For the description of the transition 
in the (/Xb-T) phase diagram, we found that the order of [ih is set for large lepton asymmetries 
by /. Assuming all particles massless and T S> Tqcd we find for lepton asymmetries 3lf = I the 

approximate solution fxb = (^^b — ij The numerical results for ^5 including all Standard 

Model particles at all temperatures with their physical masses are shown in figure [T| We find 
that sign and size of the lepton asymmetry plays an important role for the description of the 
cosmic QCD transition in the (//f,-T)-plane. For large asymmetries / < 10^^ the trajectories 
are pushed in the vicinity of a possible first order transition, predicted by results from lattice 
simulations and perturbatve calculation of QCD combined with results from heavy ion collisions 
|10j . It is a good approximation to describe the early universe and even the phase transitions 
as equilibrium processes jlT], but for a more concrete statement at T = Tqcd we would have to 
incorporate the strong interaction. To describe the cosmic QCD transition in a universe with 
baryon, lepton flavour and charge conservation, a two dimensional (/if,"T) phase diagram is not 




Figure 1. Dependence of the /i^-T trajectory on sign and size of a lepton asymmetry 



sufficient. More realistic is a five dimensional phase diagram (/u^^, , fXh, T) for a charge neutral 
and homogeneous universe. For a more detailed discussion see |6J . 

Large chemical potentials influence also the energy density and therefore the Hubble rate 
H = (8/3)7rGAr y/e{T\JjLi), with Newtons constant Gn and the energy density of all particle 
species (particle and antiparticle) e = Ej=j,i 5j7(2vr2) f^^ E^{E^ - m'j)^/^ fj{E)dE. The 

distribution function fi{E) = (exp ±1)"^ describes the boson or fermion distribution for a 
particle i. 

This influences another major event of cosmology: The freeze out of the best motivated dark 
matter candidate, the WIMP [12] with a typical mass of 10 GeV < < 1 TeV. It decouples 
statistically from the Standard Model particle^ at Tfo = 171^/25, when its annihilation rate 
drops below the Hubble rate. The number density remains constant (^x'^~^)today ~ ("'X'*~^)fo 
and thus the relative relic abundance OwiMph^ = m^n^p~^ is observable today. The number 
density follows the Boltzman equation n'^ + 3Hn^ = — {cr\v\) (n^ — n^^ gq) with the thermal 
average of the cross section, which can be approximated as (av) = a + b/x + 0{x^'^). The 
numbers a and b describe s- and p-wave annihilation and x = m^/T, for more details see 
e.g. |12j . Increasing the energy density via large chemical potentials changes the Hubble rate 
and thus the temperature, where the freeze out happens. We illustrate this effect by the change 
of the helicity degrees of freedom in the energy density (7*(T, {fii}) = ^;J^e(T, {/ii}) and in the 

entropy s{T) = ^T^gs* = - J [/In/ =F (1 ± /)ln(l ± /)] where upper and lower signs 

refer to bosons and fermions, respectively. The present relative relic abundance becomes then 
^^WIMP^^ — w2 " ( ) r,j!^h?^ ■ As long as the deviation Ag* is much smaller than the 

t^eV \ 9S* / fo <^~rO/XfQ 

standard g^{l = If = 0) we can approximate 

i'wiMP 2 V xioJ g* gs* 

with i = 0,1 depending on the s- or p-wave annihilation channel domination. For the numerical 
analysis we assume the governing effect in helicity degrees of freedom and neglect the effect in 
Xfo oc c In . This leads to: 

^wiMp{l, If / 0) 
^wiMp{l, If = 0) 

A numerical analysis with all physical particle masses and not restricted to small chemical 
potentials is shown in figure [2] for different lepton asymmetries If = We show the ratio 
^wiMp(^ ^/)/f^wiMp(^ = // = 0) as a function of lepton flavour asymmetry If and observe an 
effect of order 1 % for // = 0.01 and of almost 20% for // = 0.1. 

3. Conclusion 

In this article we have shown that large lepton asymmetries do have an impact on the evolution 
of the early universe. We focused on the cosmic QCD transition and the WIMP freeze out. 

We have shown that b <^ \l\ < 0.01 can significantly influence the dynamics of the QCD phase 
transition and maybe trigger even the order of the transition in the {fib — T)-plane. However, a 
sufficient description of the cosmic QCD transition needs at least 5 dimensions (//, 6,T, q = 0) 
and even more for an inhomogeneous universe. 

We have further shown that large lepton flavour asymmetries can decrease the relative relic 
abundance of WIMP dark matter for a given cross section. For symmetric flavour asymmetries 



l g*{i,if^o) gs*{iJf = o) 
g.{l,lf = 0) gs*{l,lf^O) 



(6) 



^ Also called chemical decoupling. The WIMPs remain in kinetic equilibrium with the Standard Model particles 
until T ~ ©(few MeV). 




Figure 2. Comparison between the relative WIMP dark matter abundance without lepton 
flavour asymmetries Vt£,M{l = ^/ = 0) s-^icl large asymmetries If = le = = It- We observe 
for If = 0.1 an effect of 20 percent. The full numerical solutions include all Standard Model 
particles. 



h = = It and If > 0(0.01) we found a few percent effect on the relic WIMP abundance. This 
might allow for smaller cross sections for a WIMP candidate to match the observed dark matter 
abundance. 
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